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Self-consistent new renormalization group flow equations for an 0(Af)-symmetric 
, scalar theory are approximated in next-to-leading order of the derivative expansion. 

^ , The Wilson- Fisher fixed point in three dimensions is analyzed in detail and various 

O ' critical exponents are calculated. 

)^ 

0^ \ 1 Introduction and conclusion 

(N 

In this talk we present novel nonperturbative flow equations for an 0(iV)-symnietric 
scalar theory in next-to-leading order of the derivative expansion. Based on the 
Wilsonian Renormalization Group method also known as the Exact Renormaliza- 
. tion Group (ERG) approach we combine a Schwinger proper time regularization 

CO \ scheme with a perturbative treatment in order to obtain self-consistent RG flow 

equations. We investigate in detail the nontrivial Wilson-Fisher fixed point in 
three dimensions and determine universal quantities such as the critical exponents 
. r], V, f3 and S independently for arbitrary numbers of field components N without 

resorting to any polynomial truncations. For this approach the typical and inherent 
analytical threshold functions in the flow equations, which describe the continuous 
decoupling of the massive modes from the RG evolution enable a direct and smooth 
Oh! linkage between the four-dimensional system at zero temperature and the critical 

0{N) universal behavior near the critical temperature Tj,. The dimensional re- 
ductimi phenomenon is embedded in the threshold functions in a very transparent 
way.ElH The prediction of the universal critical exponents is one of the great achieve- 
. ments of the ERG methods but they may also be applied to other systems e.g. with 

H ' a first-order phase transition.Q All ERG equations are equivalent up to the choice of 

the regulator or cutoff function. For instance Polchinski's equations in the limit of 
a sharp momentum cutpff are the Wegner-Houghton equations and are equivalent 
to Wilson's equations B In order to extract physical relevant quantities from the 
ERG equations which are complicated highly non-linear functional equations one 
has to employ an appropriate approximation or truncation scheme. The ERG equa- 
tions themselves are by construction scheme independent but this need not be true 
any longer in specific truncations. Thus, different realizations of the cutoff func- 
tion or regulator for the ERG equations yield an uncertainty in the determination 
of universal quantities. To leading order the critical exponents are scheme inpe- 
dendent but in next-to-leading order the results become scheme dependent. This 
makes it particularly difficult to determine uniquely the critical exponent rj. For 
instance in four dimensions the sharp cutoff Wegner-Houghton equation predicts 
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the proper field anomalous dimension 77 in next-to-leading order in the derivative 
expansion but it fails at the Wilson-Fisher fixed point in three dimensions.cl We 
propose a smooth cutoff version of the ERG based on a heat-kernel regularization 
which should overcome these difficulties and obtain remarkable precise values for 
the anomalous dimension 77 in a very simple uniform wavefunction renormalization 
approximation by neglecting the field dependence in Z^. In order to provide an 
estimate of the scheme dependence we determine several critical exponents for dif- 
ferent implementations of the smooth cutoff functions and find a rapid convergence 
to stable values that are in good agreement with entirely different approaches. An- 
other criterion for estimating the scheme dependence may be based on the minimum 
sensitivity principle which we do not further pursue in this workfl The next section 
is devoted to a brief review of the derivation of the flow equations for the potential 
Uk and the coefficient Zk of the lowest-order derivative termfl 

2 The flow equations with a proper time regularization 

As an illustrative example and as a test we apply the method with a proper time 
regulator to an 0{N) scalar theory. In the ultraviolet (UV) region of the theory at 
the scale A of the order of 1 GeV we consider the following effective Lagrangian 

u = \(d,$)\\[$^-^if (1) 

with the A^-component vector $ — ($1, $2, • ■ ■ , *&Ar)- The field $q denotes the 
minimum of the potential. The expansion in powers of momenta of the effective 
action up to order 0(9^) reads 




The perturbative one-loop contribution to the effective action yields formally a 
non-local logarithm which we regularize by a proper-time regularization resulting 
in a finite local action. We implement a multiplicative, a priori unknown, blocking 
or smearing function /^(t) which governs the coarse-graining in the proper-time 
integrand. In this way the sharp proper-time cutoff is replaced by a smooth one 
and the scale k acts as an infrared cutoff separating the low- and high-momentum 
modes of the fields. After introducing a complete set of plane wave states m the 
heat kernel this procedure yields the following d-dimensional effective actiontl 

r[$] = -1| d'^xj^ ^f,J ^tr ,~ri^~^^p.e^^e^^v:;i^)) 

with the shorthand notation V--{^) — g^J^. for the (A^ x A^)-matrix valued sec- 
ond derivative of the 0(A^)-symmetric potential V at the UV scale A. The sec- 
ond derivative is given by = A - Sij + 2X^i^j. The trace in Eq. (||) 
runs over the A^ fields and can be evaluated analytically with standard tech- 
niques. In order to find the next-to-leading order of the effective action we expand 
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exp [— T {—2ip^d^ — (3^ + V^")] in powers of derivatives up to second order. To this 
order the action separates into two parts 



r(2) 



(4) 



with r*^"^ denoting the one-loop potential contribution (containing no derivatives) 
and r^2^ the second-order contribution (containing two derivatives). Comparing the 
expansion coefficients U, Zi of the effective action (||) with the corresponding 
terms in V2 we can extract the effective potential contribution and the wavefunction 
renormalization contributions. In order to simplify the work drastically we perform 
a uniform wavefunction renormalization by neglecting the field dependence and 
consider only one wavefunction renormalization = Z^. It turns out that this 
simple approximation already describes all qualitative features at the phase tran- 
sition. The desired self-consistent flow equations are obtained by differentiation 
of the equations with respect to the infrared scale fc followed by the substitution 
V Uk on the right-hand side. This replacement has an analogy to the Schwinger- 
Dyson self-consistent resummation. The structure of the flow equations depends 
on the choice of the blocking functions but the universal physical results which are 
obtained by solving these flow equations towards the infrared should not depend 
on the specific choice. The regulator essentially governs the way how the irrelevant 
operators of the theory are integrated over. Here we chose for any integer M > 
the following form for the derivative of the blocking functions 



k 



dk 



-2(rZfcfc2)(^^+2)e- 



Zkk' 



(5) 



and obtain e.g. for the blocking function with M = 1 the rescaled dimensionless 
flow equations in d dimensions 



— 'dk — ^ ^ + (d - 2 + 77)0 + — 



I 



1 + 2u',. 



N - 1 
1 + 2u' 



(6) 



k dZk 
'Z^~dk 



2Sd 



(1 



1 



- 1 



4>=<i>a 



(7) 



where = 2/r(d/2)(47r)''/2, c^?- 



-(ci-2+r,)^2 



and Ufe(02) k-^Uk^'^^) 



and we have introduced the anomalous dimension r\ = —kdk ^nZk- A prime on the 
potential denotes differentiation with respect to cf)^. Because we have omitted the 
field dependence of the wavefunction renormalization the flow equation (^ for Z^ 
has to be evaluated at the minimum of the potential 



3 Results 



In order to test the approximations described above we solve the coupled system of 
flow equations by discretization of the field (f)"^ for a general, not truncated potential 
Uk^cj)^) numerically on a grid by a fifth-order Runge Kutta algorithm for each grid 
point. Details can be found in the work by O. Bohr et al.EI. The flow equations 
(^ and (0) are in a scale-independent form and, due to the dimensional reduction 
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Figure 1. The evolution of the dimensionless minimum towards the scaUng solution. 



phenomenon, we do not need to use the finite temperature version of the flow equa- 
tions at the critical temperature. We can directly link the 0{N) universal behavior 
near Tc with the physics at zero temperature and therefore it is sufficient to em- 
ploy only the three-dimensional zero-temperature equations in order to investigate 
the critical regime of the phase transition. In this way we circumvent the evalua- 
tion of Matsubara sums for the investigation of the critical point. A second-order 
phase transition involves an infrared fixed point of the RG transformations. Thus, 
the physics close to the phase transition is scale invariant and the critical behav- 
ior should be described by a scale-independent solution. This is demonstrated in 
Fig. |l| where the fc-evolution towards zero of the dimensionless minimum (jjQ of the 
potential for different initial values near the critical value at the UV scale k = A 
is shown. The dashed line indicates the fc-independent constant (fixed point) scal- 
ing solution while the evolution near the critical value deviates either towards the 
spontaneously broken (0g ^ 0) or the symmetric {(j)Q = 0) phase. Of course, not 
only the minimum <pQ but all quantities show a scaling behaviour. For dimensions 




Figure 2. /3a (A) for the 0(4)-model compared to a one- and three-loop e-expansion. The arrows 
indicate the flow w.r.t. the scale k towards the infrared. 
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less than four two fixed points are present in the system and can be identified by 
zeros in the beta functions. In general it is difficult to calculate the zeros since 
this requires knowledge of the physics beyond perturbation theory. In Fig. || the 
curve, denoted by RG, shows the beta function Px of the quartic coupling for an 
0(4)-model within our approach. This result is compared with a one- and three- 
loop e-expansion. Thus, the nontrivial infrared-stable Wilson-Fisher fixed point 
can clearly by identified around A* ~ 2.3 for A'' = 4. Note, that this value depends 

Table 1. Critical exponents for various TV. Our results (RG) are compared with lattice simulations, 
e-expansion (e) and a perturbativc calculation (SPS). 





V 


5 


V 


/3 


RG 

N=0 

e 

SPS 


0.039 
0.041(25) 
0.032(3) 
0.027(4) 


4.77 
4.77(14) 
4.81(2) 
4.84(3) 


0.59 
0.592(3) 
0.588(2) 
0.588(2) 


0.30 
0.308(6) 
0.304(2) 
0.302(2) 


RG 
^ lattice 
e 

SPS 


0.044 
0.044(31) 
0.038(3) 
0.032(3) 


4.75 
4.75(17) 
4.78(2) 
4.81(2) 


0.64 
0.631(2) 
0.631(2) 
0.630(2) 


0.32 
0.329(9) 
0.327(2) 
0.325(2) 


RG 
N=4 lattice 
e 


0.037 
0.0254(38) 
0.03(1) 


4.79 
4.851(22) 
4.82(6) 


0.78 
0.7479(90) 
0.73(2) 


0.40 
0.3836(46) 
0.38(1) 


N=100 RG 


0.0025 


4.99 


0.99 


0.49 


large- A' 





5 


1.0 


0.5 



on the particular form of the flow equations encoded in the choice of the blocking 
function. However, an expansion e.g. of the beta function fix in powers of the cou- 
pling yields in three dimensions /3a = —A + ^^l^^y^A'^ -f 0{}?^ which agrees with 
the first two terms of an e-expansion as it should. Corresponding to the large- iV 
limit the fixed point A* tends to zero merging with the trivial Gaussian fixed point. 
Linearizing the flow about the Wilson-Fisher fixed point allows the determination 
of the critical exponents by means of finding the corresponding eigenvalues. We 
have calculated the exponents /3, (5, v and r\ independently, which allows to test 
the well-known scaling relations among the critical exponents. Their values are 
listed in Tab. |]for various TV. The agreement with different approaches is a further 
confirmation of our approximations described above. One realizes explicitly the 
convergence of all calculated critical exponents to the large- iV values quoted in the 
last line of Tab. |^. In the large- iV limit the local potential approximation becomes 
exact and the anomalous dimension 77 vanishes. 

In order to investigate the scheme dependence of the non-truncated set of flow 
equations for the full potential we derive flow equations for the blocking functions 
defined in Eq. (||) with M = 0, 1, 2 and neglect the wavefunction renormalization. 
As an example we take iV = 4 and calculate again the critical exponent v and /3. 
The dependence of these critical exponents on the blocking function of order M is 
shown in Tab. ^ A small systematic decrease in the values with the order M is 
observed. This behaviour is also seen in the results of Refs. B'B where the dependence 
of the parameter M practically vanishes for bigger values of M. Already for M >\ 
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Table 2. The critical exponents v and fi for N = i with different blocking functions {t) 
(without wavefunction renormalization). 



M 


v 


P 





0.853 


0.42 


1 


0.815 


0.405 


2 


0.814 


0.403 



we do not see any difference in the second significant digit. Going beyond the Local 
Potential Approximation (LPA) and taking the wavefunction renormalization into 
account this tendency becomes more stable with the order M . An analysis of the 
blocking functions Z^^'* with increasing M shows ±hat the regulator has the effect 
of selecting smaller and smaller momentum shells.u It represents for very large M- 
values a kind of sharp cutoff limit which, however, cannot be directly linked to the 
sharp cutoff Wegner-Houghton equation. This means that a deeper understanding 
of the relation between this novel approach concerning the heat kernel cutoff and 
other formulations of the ERG is still an important outstanding issue. So far we and 
other authocaJiave to restrict the comparisons of the different methods to numerical 
results only.BM 
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